Figure 1: Arithmetic means of crop speculative coverage.
coverage.lm1 <- lm(speculative.pct.coverage ~ as.factor(block) + species , data = coverage, subset = (treatment != "N"))
resid_panel(coverage.lm1)
joint_tests(coverage.lm1)
## model term df1 df2 F.ratio p.value
## block 3 30 8.800 0.0002
## species 10 30 12.567 <.0001
plot(emmeans(coverage.lm1 , "species"), comparisons = TRUE) +
# scale_y_discrete(limits = rev(levels(species))) +
xlab("Speculative coverage (%)" ) +
ylab("Cover crop") +
# coord_flip() +
theme(text = element_text(size=20)) +
theme_classic()
Figure 2: Estimate speculative percent coverage of different cover crop species